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Abstract 

We find a formula for the number of permutations of [n] that have exactly s runs 
up and down. The formula is at once terminating, asymptotic, and exact. 

1 Introduction 

We will say that a run of a permutation cr is a maximal interval of consecutive arguments 
of cr on which the values of a are monotonic. If the values of a increase on the interval then 
we speak of a run up, else a run down. Throughout this paper we will use the unqualified 
term run to mean either a run up or a run down. These runs have been called sequences 
by some other authors, and have been called alternating runs by others. For example, the 
permutation 

(723851469) 

has four runs, viz. 72, 238, 851, 1469. We let P{n, s) denote the number of permutations of 
n letters that have exactly s runs. Here are the first few values of P{n, s): 

n\k 1 2 3 4 

2 2 

3 2 4 

4 2 12 10 

5 2 28 58 32 
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There is a large literature devoted to this P{n,s), which we will survey briefly in section 
EJ But although a number of recurrences and generating functions etc. are known, it does 
not seem to have been noticed that an interesting exact formula of the kind we present in 
this paper exists. Carlitz [6j has derived an exact formula for P{n, s), but that one is not at 
the same time an asymptotic formula. The Carlitz formula is discussed further in the final 
section. 

Our approach to this problem differs from previous studies in that we concentrate on the 
generating functions Us{x), defined for each fixed s > 1 by 

n>2 

whereas most earlier work has dealt with generating functions for fixed n. By finding the 
form of these generating functions we will be able to exhibit a formula for P{n, s) which is 
simultaneously 

• exact, and 

• terminating, and 

• asymptotic, for fixed s and n —>■ oo. 

To our knowledge, the asymptotic behavior of the P{n, s) has not been previously explored. 
The formula that we will find is of the form 

Pin, ^) = ^ - -^;zr- + Mn, s){s-2r + --. + z^,_i(n, s), (n > 2) (1) 
in which each ipiln, s) is a polynomial in n whose degree in n is \i/2\. 

2 Outline of this paper 

In section El we will find the generating functions Us{x) = X]n-^('^' ^)x"', as a rational func- 
tion. Since the denominator will appear in completely factored form, we can write out, in 
section m a formula for P{n, s) of the type described above. 

Interestingly, the formula will be, in that section, uniquely determined except for the 
coefficient of the leading term! That is, we will show in that section, that for fixed s we have 
P{n, s) = K{s)s'^ + . . . , but K{s) will be, for the moment, unknown. 

In sectionElwe begin the task of determining the multiplicative factor K{s). Surprisingly, 
although the tools that will have been used up to that point will be entirely analytical in 
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nature, the determination of K{s) will be done by an "almost-bijection." We will show that 
P(n, s) is, for fixed s, asymptotic to the number of s-tuples of pairwise-disjoint subsets of 
[n], each of cardinality > 2, and the asymptotic behavior of the latter is easily found. 

The combination of the former analytical results and the latter bijective argument results 
in the complete formula for P{n, s). 

3 Finding the Us{x) functions 

The recurrence formula for the numbers P{n, s) is well known and is due to Andre PP, 

P{n,s)=sP{n-l,s) + 2P{n-l,s-l) + {n-s)P{n-l,s-2), (n > 3) (2) 

with P(2, s) = 2Ss,i- From this recurrence one finds easily a recurrence for the generating 
functions Us{x) = ^^P(n, viz. 

(1 — sx)us{x) = 2xUs^i{x) + x^u'g_2{x) — (s — l)xus^2{.x)i > 2) (3) 

with Ui{x) = 2x^/(1 — x), Uq{x) = 0. The next three of these functions are 



U2{x) 
■U4(x) 



(l-x)(l-2a;)' 

2x^(5 - 6x) 
(l-3x)(l-2x)(l-x)2' 

4x^(8 - 29x + 24x2) 
(1 - 4x)(l - 3x)(l - 2x)2(l - x)2' 



We will find the general form of these functions, and from that will follow the desired 
formulas for P{n, s). 

Theorem 1 We have, for each s = 1, 2, 3, . . . , 

""'^"^^ " (1 - sx)(l - (s - l)x)(l - (s - 2)x)2(l - (s - 3)x)2 ■ • ■ (1 - x)L(«+i)/2J ' ^ ' 



where $s(x) is a polynomial of degree 1 + 



s(g+2) 



The degree of the denominator is 



4 

which is exactly 1 less than the degree of the numerator, for all s > 1. 



s{s+2) 
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3.1 Proof of Theorem [T] 

The proof is by a straightforward, though tedious, substitution of the form Q into the 
recurrence (jHl) to find a recurrence for the numerator polynomials This will establish 

that they are indeed polynomials and will provide the claimed degree estimates. We will do 
this by putting every term over the common denominator 

A,{x) = (l-sx)(l-(s-l)x)(l-(s-2)a;)2(l-(s-3)x)2---(l-a;)L(^+^)/2J 

s-l 

= 11(1 -(^-^)^)"' 

?:=o 

where we have written {ei}j>o = {1, 1, 2, 2, 3, 3, 4, 4, . . . }. 

For technical reasons it will be useful to rewrite the recurrence (0) in the form 

, , 2xUo_i(x) x'^u' r,(x) (s — l)xUs-2ix) 

u,{x) = , ^ ' + , - ^ — — 

(1 — sx) (1 — sx) (1 — sx) 

(1 - sx)As_i(x) (l-sx)As_2 (1 - sx)As_2 As_2(a;) 

(g - \)X^S-2{X) 
(1 - SX)[\s-2{x) 

1 \ 2x^,.x{x)A,{x) x'^',^^{x)As{x) x'^,^2{x)As{x)^,_^{x) 
A,(x)l(l-sx)A,_i(x) (l-sx)A,_2 (l-sx)A,_2 A,„2(a;) 
(g- l)x'l',_2(x)A,(x) ) 

(1 - gx)A,_2(x) j ^ ^ 

Each of the four terms inside the braces is a polynomial in x whose degree we will now find. 
Consider the ratio 

A,(a:) ^ n;:;(i-(^-j>^' 

(1 - gx)A,_i(x) (1 - sx) n;:o(l - (s - 1 - 

(i-^^)n;:;(i-G^-jw-^ 

j=l J even; 2<j<s— 1 

which is a polynomial of degree [(s — 1)/2J. 
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It follows that 



As(x) 



(1 - sx)As-2{x) 



AJx) 



As-i{x) 



(l-(s-l)x) 



(1 - sx)As-iix) J \{l-{s- l)x)A,_2(x), 

J even; 2<j <s— 1 j even; 0<j<s— 2 

n n a-is-j)x) 

J even; 2<j<s— 1 iodd;l<j<s— 1 

is a polynomial in x of degree s — 1. 

We can now deal with the third term of the four inside the braces in . Since 



A' 



s-2 



X 



A,-2ix) 



s-1 



i=2 



ej-2{s-j) 



- (s - J>' 



we have 



x^^s-2{x)A,{x) Al 



s-2 



X 



(1 - Sx)As-2 As-2{x) 



X^(l>s^2ix) 



's-1 

Y[{l-{s-j)x) 



' s-1 



0=2 



ej-2(g - j) 
- (s - j)x 



If d{s) denotes the degree of then this last member is a polynomial in x of degree 

2 + d{s - 2) + s - 2 = d{s - 2) + s. 

We have now shown that each of the four terms inside the braces in (0) is a polynomial 
in X. Their respective degrees are 

d(s - 1) + 1 + [{s - 1)/2J , d{s -2) + s, d{s -2) + s, d{s -2) + s. 

Hence we have d{s) = max {d{s — 1) + [{s + 1)/2J , d{s — 2) + s), with d{2) = 3 and d{3) = 5. 
It is remarkable that this difference equation has a simple solution. Its solution is 



d{s) 



1 + 



s{s + 2y 



as can easily be checked, and in fact all four terms inside the braces in (jH)) have the same 
degree! This completes the proof of the Theorem. 
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4 The formula for s) 

From the partial fraction expansion of ^ we find at once that 

P{n, s) = Mri, + ^i(n, s){s - 1)" + Mn, - 2)" + ■ ■ ■ + s), {n > 2) (6) 

where each '?/'j(n, s) is a polynomial in n of degree at most [V2j , and it remains to find these 
polynomials. To do this we substitute © into the recurrence ^ and match the coefficients 
of each term (s — i)". The result of this substitution is that 

(s - i)ipi{n, s) = sipi^n - 1, s) + 2^/'j_i(n - 1, s - 1) + (n - s)^/'i_2(n - 1, s - 2). 

Perhaps the best way to find these ip^s explicitly is to assume a solution in the form of a 
polynomial in n of degree and solve for the coefficients of that polynomial. We can 

begin with tp-i{n, s) = and ipoi^, s) = K{s) (since ipo is of degree zero in n) where K is to 
be determined. We then find that 

Mn, s) = -2K{s - 1), 7/'2(n, s) = ^K{s - 2)(s + 8 - 2n) 

ijsin, s) = Ik{s - 3)(2n - s - 3), M^, s) = ^K{s - 4){4n^ - 4n(s + 8) + + 15s + 32) 
For example, for s = 4 we find 

1/4 (-1) (-1/2) 7/2 2 (-9) 

u,{x) = + + ; , + -\- + ^ + + 2x + 19/4. 

1 — 4x 1 — 3a; (1 — 2xy 1 — 2x [1 — xy 1 — x 

From this it follows 

P(n, 4) = 4"-^ - 3" + (6 - n)2"~i + {2n - 7), (n > 2). 



5 The factor K{s) 

We have now described the formula for P(n, s) completely except for the multiplicative factor 
K{s). It remains to show that K{s) = 2~^^~'^\ For this, it would suffice to prove the next 
Theorem for fixed s and n ^ oo] since the proof is applicable to a larger range of s, we state 
it in that manner: 

Theorem 2 Let e > 0, and {(n, s)} he an infinite sequence of pairs such that n —>■ oo and 
s < (1 + e)^"*^?!/ logn. Then, 

Pin,s)r^^^s\ (7) 
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5.1 Proof of Theorem [2] 

To fix ideas, we will do this by showing that the number P{n, s) of permutations of n letters, 
with s runs, the first of which is a run up, is ~ s"/2''^^. Evidently the number for which 
the first run is down will be the same, and the desired result will follow. Henceforth we will 
always assume that the first run is a run up. There are two steps to the proof. In the first 
step, we show that the set of permutations counted by -P(n, s) can be put into bijection with 
s-tuples of subsets {Si, . . . , Sg) (each Si C [n]) satisfying certain properties. In the second 
part of the proof, we introduce a function called $ whose domain is the Cartesian product 
of these s-tuples with a set of cardinality 2*~^, and whose range is a set of size s". We prove 
that this function $ is an injection. Although we have no succinct description of the image 
of this injection, we are able to show that for (n, s) in the range hypothesized by the theorem 
the image is asymptotically all of the range set. 



5.2 First part of the proof 

Let n(n, s) be the set of all n-permutations with s runs up and down, the first of which is 
up. Let n(n, s) be the collection of all s-tuples (^i, . . . , Ss) of nonempty subsets of [n] which 
are almost pairwise disjoint, in that 



\Si^Sj\ 




Further we require that 
and that 



i + 1 and 1 < i < s; 



\S^\>2, Vz, (9) 



max {Si) = max (S'j+i) E SiH S'j+i (V odd i) 

min (Si) = min (Sj+i) G Si fl Si^i (V even i). (10) 



Lemma 1 The number of s-tuples of subsets of [n] that satisfy /f^)-/ li(^) is equal to the 
number of permutations of [n] with s runs, the first of which is up. 

Indeed to reconstruct the permutation from the s-tuple of sets, we first sort each of the sets, 
the first in increasing order, the second decreasing, etc., then merge the sets, and finally 
delete one element of each of the adjacent duplicates that appear. □ 

Hence it suffices to show that the number of s-tuples of subsets that satisfy (jHjl- ()lUp is 
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5.3 Defining the function $ 

By a choice sequence h = {hi, . . . , hg-i) we mean an s — 1-tuple where each hi is either 
equal to z or to i + 1. The set of all such choice sequences will be Hg. The function to be 
constructed is a mapping 

$ : /f, X n(n, s) {(Ti, T2, . . . , T,) : Vz, T, C [n]}. 

Let h G i^s, and let (5*1, . . . , S'^) be a family of subsets satisfying (Pl- lfTUj) . For each 
i = 1, . . . , s — 1, let Ci be the unique element that belongs to Si Si^i. These Cj's are all 
different, since = Cj with i < j would imply that Si fl 5^+1 is nonempty, contradicting 
(jSJ. Perform the following s — 1 delete operations: for each i = l,...,s — 1, delete the 
element Cj from the set Sh^- The resulting s-tuple of sets remaining after these deletions is, 
by definition, $ (h, {Si, . . . , Ss)) ■ 

The image of this mapping does not include all s-tuples of sets, as the following Lemma 
shows. 

Lemma 2 // (Ti, . . . , Tg) is in the image of $ then 

1. the Ti's are pairwise disjoint, and 

2. the union of the Ti 's is [n] . 

□ 

It is possible for some of the Tj's to be empty. We remark that the number of s-tuples 
(Ti, . . . , Tg) in which the Tj's are pairwise disjoint and whose union is [n] is s^. 

5.4 The mapping $ is injective 

The way we prove this assertion is to give a reconstruction algorithm. The algorithm begins 
with an s-tuple (Ti, . . . , Tg) of subsets which putatively belongs to the image of $. It attempts 
to reconstruct the preimage. It will be clear from the algorithm that the preimage can be 
only one thing, if it exists at all. There is one "early exit" point in the algorithm where the 
search for a preimage is abandoned, because it obviously does not exist. If the algorithm 
executes all the way to finish, then it will have found the only possible candidate for a 
preimage. However, it is still possible that the s-tuple of sets found at the end will not 
satisfy one of the required conditions (jH|)- (piH) . 

Lemma 3 The mapping $ is injective. 
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Proof. Let (Ti, . . . ,Ts) be an s-tuple of pairwise disjoint (possibly empty) sets whose union 
is [n]. Here is the reconstruction algorithm: 

1. (Find consecutive unions) It is easy to see that if the T's are in the image of $, then 

Ti U Tj+i = SiU Si+i, 1 <i < s. 
So if one of the inequalities 

|TinTi+i|>3 l<i<s, 

fails, then the reconstruction fails and no preimage exists. Otherwise we can reconstruct 
all of the unions Si U Si+i. 

2. (Reconstruct the set of deleted elements) Put ei = max {Si U S2), 62 = min (5*2 U 5*3), . . . . 

3. (Recover the choice sequence h) For each i = 1, . . . , s — 1, since Cj G Tj U Tj+i, and 
because the Tj's are pairwise disjoint, there will be exactly one index, hi, say, such that 
hi e {i,i + 1} and ^ T^,. 

4. (Re-insert the elements that were deleted) For each i, 1 < i < s, insert the element Cj 
into the set T/^.. 

If the reconstructed sets {Si, . . . , Sg) satisfy (injl-lfTUj) then we have found the unique preimage. 
Otherwise no preimage exists. □ 

Thus if P{n, s) is the number of permutations of n letters with s runs, the first of which 
is up, then we have shown that 

2'-^P{n,s) < s"". (11) 

5.5 When does the algorithm terminate without a preimage? 

If the reconstruction algorithm does not early exit in step 1, yet fails to find a preimage, 
then one of the conditions (jSjl- lfTUI) is not satisfied. We will now visit each of these in turn 
to see when it might fail 

1. (Can (jHI) fail?) The intersections TjflTj+i were all empty before the insertions; however, 
the operation, "insert element Cj into Th" either added an element of Tj to the set Tj+i, 
or vice- versa. That operation alone caused the two adjacent sets to have intersection 1. 
The only other insertion which could have affected Tj is the one which involves element 
ej_i. If that operation increased the size of Tj, then it did so by inserting an element 
from Tj_i, which element could not possibly be present in Tj+i. Thus, the only other 



9 



insertion which could possibly affect the set Tj will have no effect on the cardinality of 
Tj n Tj+i. Likewise, the only other operation which can possibly affect the cardinality 
of Tj+i will have no effect on the cardinality of Tj fl Tj+i. So, the intersection Si fl Sj+i 
will always have size 1, as required. 

2. (Will Si n S'j = when j > i + 1?) Only the case j = i + 2 is not obvious. If Si fl S'j+2 
is not empty, then during reconstruction some element originally belonging to Tj+i 
was inserted into both T, and Tj+2. (Any element originally in T, cannot end up in 
5*4+2, and vice- versa.) This means that some element e G Tj+i is both the maximum of 
Tj U Tj+i, as well as the minimum of Tj+i U Tj+2; (or the other way around). But 

max (Tj+i) < max (Tj U Tj+i) = Cj = min (Tj+i U Tj+2) < min (Tj+i), 

so Tj+i has just one element, Cj, and that element lies between the surrounding sets Tj 
and Tj_|_2- This is how the reconstructed sets can fail to satisfy (jH)). 

3. (Can Q fail?) Yes, if one of the sets Tj has cardinality or 1, then it is possible that 
the cardinality of the reconstructed Si may be less than 2. 

4. (Can (fTUI) fail?) No. By the nature of the reconstruction, the S'j's always have this 
property. 

We can now prove 

Lemma 4 // in the given sequence T = (Ti, . . . ,Ts), all sets have cardinalities at least 2, 
then T has a preimage under $. 

For then the unions TjUTj+i have size 4 or more, so we don't terminate the reconstruction at 
Step 1. The only other two possible failures — when an intersection Sj nS'j+2 was nonempty, 
or one of the Si was too small — could both be traced back to a set Tj which had size or 
1. □ 

A crude lower estimate from Bonferroni's inequalities tells us that the number of s-tuples 
T that are pairwise disjoint, with union equal to [n], and with all cardinalities > 2 is at least 

s"-(n + s)(s-l)"-\ 

The reason: s{s — 1)" is an upper bound on the number of T's for which some component 
is the empty set; and ns{s — 1)"^^ is an upper bound on the number of T's for which some 
component has cardinality one; then, s{s — 1) + ns < s{n + s). Hence 

2'-^P{n, s) > - s{n + s){s - l)"-\ (12) 

which, taken together with (lllj) completes the proof of (|7j), since our hypothesis on the pairs 
(n, s) implies 

(n + s)(s - 1)"-^ < 2n(s - 1)""^ = o(s"). 
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6 Survey of the literature 



Andre was the first to study pj the runs up and down of permutations, and the fundamental 
recurrence (j21) is due to him. His paper includes a table of P{n, s) through n = 8, with one 
error in the final row. A great deal of information about P{n, s) is found in vol. 3 of jH] 
(see particularly ex. 15, 16 of sec. 5.1.3). Comtet [8^ devotes an extended exercise, see page 
260, to the topic. The two variable generating function given there, however, is incorrect. A 
correct version appears in the discussion accompanying sequence A059427 of [TUj . 

Carlitz El EI visited this subject several times. In [S] he gives the two-variable gener- 
ating function 



and in jB] he finds an explicit formula for P{n, s) and information about an associated poly- 
nomial sequence. There is something wrong with the final formulas of this latter work, 
however; these formulas suggest P(8, s) = 0, 2, 250, 2516, 7060, 7562, 2770; whereas, in fact, 
P(8,s) = 2,252,2766,9576,14622,10332,2770. (Empirically, his formula always gives the 
right value for P{n,n — 1).) Further evidence that something is amiss concerns the auxil- 
iary quantity Knj] the summation formula given for this quantity does not give the values 
displayed in the table. Whether the problem can be easily repaired, we have not investigated. 

More recently. Bona and Ehrenborg |lj have proven log-concavity: P(n, s)^ > P{n,s — 
l)P{n,s + 1). In the later book [Sj, the stronger assertion, that Pn{x) = P(?2, s)x^ has 
all its roots real and negative, is made. A proof of this can be based on the relation 



which itself is a consequence of the basic recursion (j21). This implies, once it is established 
that the variance becomes infinite with n, that the numbers P(n, s) satisfy a central limit 
theorem. (That is, are asymptotically normal.) Due to log- concavity, one may deduce (see 
Theorem 4 of [2]) a local limit theorem. This leads to an asymptotic formula for P(n, s) for 
s in a different range than in our Theorem |21 
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